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ATTACK TO DETERI'dINE THE LONGITUDINAL STABILITY AND 
CONTROL CHARACTERISTICS OF THE SIDEWINDER MISSILE 
AT MACH NUMBERS FROM 1.2 to 2.1. 

By Clarence A. Brown, Jr. 


SUMMARY 


The longitudinal aerodynamic characteristics of the SIDEWINDER 
missile have been determined in flight with a full-scale rocket-powered 
model by the Langley Pilotless Aircraft Research Division. Static and 
dynamic longitudinal stability and control derivatives of this inter- 
digitated canard-wing missile configuration were determined, using the 
pulsed-control technique, at low angles of attack and for a Mach number 
range of 1.2 to 2.1, 

The lift-curve slope for the SIDEWINDER missile showed only small 
nonlinearities with changes in control deflection or angle of attack 
but indicated a difference in lift-curve slope of approximately 7 per- 
cent for the two control deflections of B ^ 5«0° and B ^ The 

large tail length of the missile tested was effective in producing 
damping in pitch throughout the Mach number range tested. The aero- 
dynamic-center location was nearly constant with Mach number for the 
two control deflections but was shown to be less stable with the larger 
control deflection. The increment of lift produced by the controls 
was small and positive throughout the Mach number range tested, whereas 
the pitching moment produced by the controls exhibited a normal trend 
of reduced effectiveness with Increasing Mach number. The effectiveness 
of the controls in producing angle of attack, lift, and pitching moment 
was good at all Mach numbers tested. 
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HraRODUCTION 


At the request of the Bureau of Ordnance ^ Department of the Havy^ 
the Langley Pilotless Aircraft Research Division has initiated a program 
to investigate the general aerodynamic characteristics of the Kaval 
Ordnance Test Station's SIDEWINDER missile. The model used in the flight 
test presented herein was a full-scale^ rocket-propelled test vehicle. 

This paper presents the results from a flight test investigation using 
the pulsed-control technique to determine the static and dynamic longi»> 
tudinal stability and control derivatives and drag data for a canard- 
missile configuration. The methods for obtaining these data are pre- 
sented in references 1 and 2. This investigation was conducted at 
a small angle-of -attack range and for a Mach number range of 1.2 to 2.1. 

The model used in this investigation was flight-tested at the Langley 
Pilotless Aircraft Research Station at Wallops Island^ Va. 


SYMBOLS 


S^ total wing area in one plane including body intercept;, sq ft 
Sg exposed canard area in one plane, sq ft 
c • wing mean aerodynamic chord, ft 


S^ body cross-sectional area, sq ft 

d body diameter, ft 

W model weight, lb 

2 

ly moment of inertia about Y-axis, slug-ft 

Ijj moment of inertia about X-axis, slug- ft? 

R Reynolds number 

M Mach number 

g acceleration due to gravity, ft/sec? 

q dynamic pressure, Ib/ft? 
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b exponential damping constant^ per second 

P period of oscillation, sec 



A aspect ratio 

a angle of attack, deg 

“brim angle of attack, deg 

a 57^ radians/sec 


6 control deflection, deg 

0 pitching velocity, radlans/sec 

An/g normal accelerometer reading, g imits 

Aj/g longitudinal accelerometer reading, g units 

A^/g transverse accelerometer reading, g units 
Cjj drag coefficient, A^/g cos a + A^/g sin 

lift coefficient, ^Aj^/g cos a + Aj/g sin 


Cm 

ACl 

AS 

AS 

Cn 

Cy 


w 

qSb 


w 

pitching-moment coefficient. Pitching moment about , center of . gravity 

SiSi^d 

average lift coefficient per unit control deflection 

average pitching-moment coefficient per unit control deflection 

normal-force coefficient, A_,/g 

qSb 


lateral-force coefficient. 




"R 


resultant-force coefficient corrected for trim. 


“ ^Ntrim)^ 
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Derivatives: 

SCt 

= _-i!^ per degree 


a 




Cjn = per degree 

cc ^ 




per degree 

6 S6 


da 


'm 


"m. 


per radian 


‘q ./0d' 

„ SCm 

P®^ radian 

a ^fcbd] 

® a? 


MODEL AND APPARATUS 
Model Description 


Sketches of the rocket-powered model used in this test are shown in 
figure le Sketches of the canard surface and wing surface are shown in 
figure 2. Photographs of the model and model booster combination are 
shown in figures 3 and 4 . Physical characteristics determined by preflight 
measurements are presented in table I. 

The body of the model had a maximum diameter of 5 inches with a 
fineness ratio of 22.95» The nose section consisted of a 2»6-lnch»radius 
spherical segment that was faired into the 5“lnch~dlameter body. Pro- 
truding in front of the spherical nose section was a sting used to mount 
part of the models instrumentation (figs. 1 and 5)« The canard surfaces 
were of arrow wing plan form with a modified hexagonal airfoil section 
having a maximum thickness at the wing -body juncture of 3-3 percent 
(fig, 2), The leading edges of the canard surfaces were swept back 66°3^’ 
and the trailing edges were swept back l4°51* (fig- 2). The canard sur- 
faces were pivoted about a hinge line located at 46 percent of the mean 
geometric chord. 

The wings were interdigitated 45° to the canard surfaces and were 
of trapezoidal plan form with the leading edge swept back 45° (fig- 2). 

The wing had a modified hexagonal airfoil section with a constant thick- 
ness corresponding to a thickness ratio of 1.2 percent at the wing -body 
juncture . 
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The control surfaces were actuated hy a slow-acting solid-propellant 
gas “driven servo explained in reference 3» In order to use the gas-driven 
servo to abruptly disturb the model in pitchy it was necessary to reduce 
the rise time of the control surfaces* Static firing of the pulse unlt^ 
prior to flight test of the model and at the desired control deflection 
range, resulted in changes that reduced the rise time to approx±mately 
0.06 second. The length of time for a fixed control deflection was 
approximately 0.6 second throughout the flight of the model. The solid 
propellant used in the gas-driven servo was capable of operating the 
pulse mechanism for approximately 10 seconds* 


Instrumentation 

The model was equipped with an HACA nine-channel telemeter which 
transmitted a continuous record of normal (two locations), transverse 
and longitudinal accelerations, angle of attack, angle of sideslip, 
control deflection, total pressure and static pressure. The transverse, 
longitudinal, and one normal accelerometer were located so as to be on 
the center of gravity of the model when the sustainer motor had burned 
outj and a second normal accelerometer was moionted on the model center 
line and inches ahead of the center of gravity of the model. Angle 
of attack and angle of sideslip were measured by a free-floating vane 
mounted on a sting which protruded from the nose of the model. Total 
pressure was obtained by a total-pressure tube extended from the fuselage 

ahead of the wings and in a plane 22 ^° to the main wing and canard 

surfaces, A static-pressure orifice was located on the cylindrical sec- 
tion of the fuselage ahead of the canard surfaces. Approximate values 
of rate of roll were obtained by MCA spinsonde equipment in conjunction 
with the telemeter antenna which was plane polarized. 

Model velocity was obtained from the CW Doppler radar and the model 
trajectory was determined through use of an MCA modified SCR-584 radar 
tracking unit. A radiosonde, released at the time of flight, was used 
to obtain atmospheric data throughout the altitude range traversed by 
the model. 


TEST TECHNIQUE 


The model was launched at approximately ^5° elevation angle from 
a zero-length launcher as shown in figure 4, The model was boosted to 
supersonic velocity by two 6~lnch— diameter solid-propellant rocket motors 
which together delivered approximately 12,000 pounds of thrust for 
3.0 seconds. After separation from the booster, a sustainer motor, made 
as an integral part of the model, delivered approximately 2,500 pounds of 
thrust for 2.6 seconds and propelled the model to the peak Mach number 
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of 2.h’J. Following the sustainer burnout^ the model was disturbed in 
pitch by a programmed square“Wave deflection of the canard surfaces « 
Transient responses to the step input of the control surface were con- 
tinuously recorded in the form of time histories as the model decelerated 
through the Mach number range « 

The canard control surfaces in the horizontal plane were pulsed in 
a square-wave motion by deflecting them abruptly to a 6 ~ +5»0° and 
holding them in that position for a predetermined time interval of 
approximately 0.6 second, then deflecting them again abruptly to a 
6 ~ -0,3° and holding them again at this deflection for the same time 
interval. This sequence was repeated throiighout the flight of the model. 
The canard control surfaces of the vertical plane were locked in the 
zero position and remained in that position throughout the flight of the 
model. 


PRECISIOM OF DATA 
Correction 


The velocity data, as obtained by the CW Doppler veloclmeter, were 
corrected for flight^path curvature and wind effects at altitude. The 
magnitudes and direction of these winds were deteimined by tracking the 
radiosonde balloon. 

In order to obtain the angle of attack at the center of gravity of 
the model, the angle of attack measured at the nose was corrected for 
model pitching velocity by the method presented in reference 4 . 


Accuracy 

The maximum possible errors in accuracy of the quantities listed 
below, on the basis of the accuracies of the instrumentation and dynamic 
pressure are presented for two Mach numbers j 


M 

Limit of accuracy of — 

M 

a 

6 

Cl 

Qm 

^^in 

1,30 

+ 0.01 

+0.50 

+ 0.20 

±0.17 

+0.60 

+0.08 

2.00 

±0.02 

±0.50 

+0.20 

±0.06 

+0.16 

+0.05 
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It should be pointed out that the tabulated quantities listed in the 
table on the preceding page are based on body cross-sectional area. From 
a consideration of previous experience^ probable errors are 50 percent less 
that those just quoted. Parameters dependent upon differences in measured 
quantities or slopes such as are much more accurately determined than 

the previously mentioned errors would Indicate, 


RESUnr MD DISCUSSION 


Complete data were received for the model tested for a Mach number 
range of 1,2 to 2,1. The Reynolds number of this test ranged from approx- 
imately 4,7 X 10^ to 10.7 X 10^, per foot, Variation of Reynolds number 
with Mach number for this test is shown in figure 5 . 


Lift Coefficient 

Shown in figure 6 are plots of lift coefficient against angle of 
attack for the two control deflections of 6 » -0,3° 6 « -!- 5»0 . 

These plots are typical lift coefficient against angle-of-attack plots 
for this model. In order to present lift coefficient against angle of 
attack and retain clarity^ only one increasing and decreasing angle of 
attack for each disturbance is presented. The hysteresis noted in the 
data is not unusual for a canard-type configuration. Several other 
pulsed control models have also experienced this hysteresis characteristic^ 
(refs, 1 ;, 2 ^ and 5 ) and the effect upon the lift-curve slope is negligible. 
Presented in figure 7 are the lift-curve slopes against Mach number from 
figure 6 and for similar plots at other Mach numbers. The lift-curve 
slope showed only small nonlinearities with changes in control deflection 
or angle of attack but indicated a difference in lift-curve slope of 
approximately 7 percent for the two control deflections of 6 + 3 » 0 ° 

and 5 « -0,3°. As would be expected;, the lift-curve slope exhibited a 
smooth variation with Mach number throughout the Mach number range 
tested. 

Also presented in figure 7 is a theoretical lift»curve slope cal- 
culated from reference 6 . Agreement between the theoretical and meas- 
ured lift-curve slopes below a Mach number of 1,60 is poor. Some of 
this disagreement might be explained {in tliat for the configuration 
tested the low aspect ratio of the rearward surfaces and reduced Mach 
numbers decreased the effective wing aspect ratio (pA) to a value less 
than one. When the effective aspect ratio is less than one the basic 
theoretical lift-curve as determined by linear theory for reference 6 
was beyond the scope of the theory and necessitated using a theoretical 
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lift-curve slope wliich the authors of reference 6 found necessary to 
extrapolate. Above a Mach number of 1.6o the agreement between the meas 
ured lift-ciirve slope and the theoretical values varied from 4 to 
15 percent « 


Dynamic Stability 

The exponential damping constant b is presented in figure 8 for 
the two control deflections . The damping-in-pitch derivative 
Cm t Cm® obtained from the faired curve of b is presented in figure 9* 
q a 

The damplng-in-pltch derivative 0™ + increased from -7^000 at a 

q a 

Mach number of 1.24 to -llj820 at a Mach number of 1.55^ then decreased 
gradually to a value of “8,600 at a Mach number of 2.15. As a result 
of the model roll rate being between 0 and 5 radians per second^ it was 
necessary to analyze the resultant~force coefficient time history of the 
normal and transverse motion by the method presented in reference 2 to 
obtain the damping of the model. 

The damping~in>»pitch derivative C™ + 0^. for the model of the 

q a 

present test compares favorably with the damping of the model of refer- 
ence 5 and^ as might be expected, the large tail length of the present 
model was effective in producing damping in pitch for the Mach number 
range tested. 

Static Stability 

The longitudinal period of oscillation of the model using the 
resultant force time histories is presented in figure 10 as a function 
of Mach number. 

Two methods were used in obtaining the pitching-moment derivative 

Cm presented in figure 11. The faired curve of Cj^ was reduced 
a a 

from the faired curve of period of oscillation of the model. The plotted 
points of figure 11 were obtained by taking the slopes of pitching 
moment against angle of attack. Plots of the pitching moment against 
angle of attack are presented in figure 12 for the two control deflec- 
tions of 6 « - 0 . 5 ° and 6 « +5.0°, The total pitching moment was 
obtained by the use of the two normal accelerometers, one located at 
the model center of gravity and the other located 4-5 inches ahead of 
the model center of gravity. The part of the pitching moment due to 
the angle of attack can then be obtained by subtracting that part which 
was contributed by the model damping. As may be seen in figure 12 some 
hysteresis is noticeable for most of the control deflections and 
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examination of these curves reveals only slight nonlineaxities . Com- 
parison of these slopes and the pitching«.moment derivative obtained from 
the period of oscillation, figure 11, indicates good agreement between 
the two methods of obtaining this derivative. 

Aeordynamio»center location was determined from the curve and 

the faired Ct curves and is presented in figure 15 in terms of inches 

from station 0 against Mach mmiber. The two curves for the aerodynamic- 
center position resulted from the two values of for the two control 

deflections. The aerodynamic center was nearly constant with Mach number 
for either control deflection but showed to be slightly less stable with 
6 ~ +3*0° than with 6 » -0.5° « 

Also included in figure 13 are the loaded and empty center ~of -gravity 
location in station niimbers and a theoretical aerodynamic-center location 
from reference 6, As mentioned previously, seme of the disagreement 
between the experimental and theoretical values is due to using a theo- 
retical llft~c\nrve slope, determined by linear theory for reference 6, 
which was beyond the scope of the theory and necessitated using a theo- 
retical lift-curve slope which the authors of reference 6 found necessary 
to extrapolate. 


Control Effectiveness 

The trim angles of attack are shown in figure ll<- as a function of 
Mach n'umber for the two control deflections of & « +5*0° and 
6 ~ “0.3°. It was expected that for the symmetrical model tested, the 
model would trim at a = 0 for 5=0. The apparent out of trim shown 
in figure l4 may have resulted from indicator out of trim or asymmetries 
due to model construction. 


The ratio of the trim angle of attack to the trim control deflection 
for the model tested is presented in figure 15. The effectiveness of 
the controls in producing angle of attack can be seen in figure 15 to 
decrease slightly with increasing Mach number until M = I.60, then 
increase gradually with Increasing Mach number. It should be noted that 
the effectiveness of the controls in producing angle of attack is as 
good at low supersonic Mach nimibers (M = 1.3) as at the higher supersonic 
Mach nijmbers (m = 2 .l), 


The control derivatives 




A5 


and 


ACm 

A6 


presented in figures 16 and I7 respectively. 


of the model tested are 
Figure 16 shows that the 


increment of lift due to the canard deflection 


A 6 


is small and 
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positive at all Mach numbers tested. The fact that is positive ^ 

means that the lift of the control surface itself exceeds the loss of 
lift on the wing due to downwash of the canards. This can prohahly he 
attributed to the interdigitation of the wings and canards reducing the 
downwash at the small angles of attack of the test and to the size of 
the canards themselves. Presented also in figure l6 is a plot of Cr 

as obtained from reference 6, Both reference 6 and the test 

A o 

indicated positive values throughout the Mach number range tested. 

Pitching effectiveness of the model (fig. I?) positive through 
the Mach number range tested and decreased from 1.^6 at a Mach number 
of 1.50 to a value of I.08 at a Mach number of 2.10 exhibiting a normal 
trend of reduced effectiveness with increasing Mach number. 

For the large static margin of the test^ 16 to 21 inches,, the effec~ 
tiveness of the controls in producing angle of attack;, lift, and pitching 
moment is good. Since this particular configuration will fly near trim 
conditions at all times, control surfaces such as these will give good 
maneuverability and still remain near trim conditions. If greater maneu” 
verability is desired, hot/ever, it is possible to reduce the static margin 
approximately 7 inches and still retain a stable configuration at low 
supersonic Mach numbers and loaded conditions. 


Drag 

Drag data for the model tested are presented in the form of Cjj ? 

min 

based on fuselage cross-sectional area, against Mach number in figure I8, 
The mlnlmimi drag coefficient varied smoothly with Mach number from l.ii-8 
at a Mach number of 1,15 "to 1.20 at a Mach number of 2 , 20 , 

Also presented in figure I8 are the zero-lift drag coefficient 
against Mach mmiber for two SIDEWINDER configurations (models of refer- 
ences 7 and 8) and the C-n for the model tested. The zero~llft drag 

nnln 

coefficient for the model in reference 8 was not published in that report 
but is presented here for comparison. The models of references J and 8 
were similar to the model tested with several changes made to the models 
that affect the drag coefficient. These changes were as follows? 

1 , Models of references J and 8 did not have an angle-of-attack, 
angle-of-sideslip indicator, 

2 , Models of references 7 and 8 were approximately 8 inches shorter, 

5. Models of references 7 and 8 were roll-rate stabilized with air- 
driven gyro-actuated rollerons. The model of reference J experienced a 
high frequency roll instability in flight which caused the rollerons to 
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deflect violently at the same frequency as the roll instability » Ihe 
model of reference 8 eliminated this roll instability without changing 
the exterior parts of the model. Although the model of reference 7 
experienced a high frequency roll instability^ the roll-rate stabili- 
zation system did roll»rate stabilize the model within a roll rate of 
+20° per second. Comparison of the drag coefficient of the model of 
reference 7 the model of reference 8 shows that elimination of the 
roll instability thereby eliminating the violent deflecting of the roUerons 
resulted in a reduction of drag coefficient of approximately 8 to ^4-0 per- 
cent at Mach numbers of 1,95 1«25^ respectively. 

As previously mentioned, the model used in this Investigation did 
not have a roll-rate stabilization system as did the models of refer- 
ences 7 8 but an angle-of -attack, angle-of-sldeslip indicator pro- 

truded from the nose of the model. The drag coefficient of the present 
model tested and the model of reference 8 differed from approximately 
20 percent to lO percent at Mach numbers of I.I 5 and 2,18, respectively. 

This reduction of drag coefficient can be attributed in part to the 
elimination of the roll-rate stabilization system and the angle-of- 
attack, angle-of-sideslip Indicator acting as a windshield or spike. 
Experimental results have shown that drag reductions to spherical nose 
sections can be made by the use of spikes or windshields. For example, 
reference 9 shows that the addition of a spike to a spherical nose will 
reduce the drag of the model from I 5 percent at a Mach number of 1.30 to 
40 percent at a Mach number of 1,60, 


COWCDJSIONS 


The results of a flight test of the SIDEWIKDER missile for a Mach 
number range of 1.23 to 2,1 indicated the following conclusions! 

1. The lift-curve slope showed only small nonlinearities with 

changes in control deflection or angle of attack but indicated a dif- 
ference in lift-curve slope of approximately 7 percent for the two con- 
trol deflections of 5 +3*0° and 6 « -0,3°. 

2. The large tall length of the model was effective in producing 
damping in pitch throughout the Mach number range tested, 

3 . The aerodynamic-center location was nearly constant with Mach 
number for either control deflection but showed to be slightly less 
stable with the larger control deflection. 

i|-. The increment of lift produced by the controls was small 

and positive throughout the Mach number range, 
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5= The effectiveness of the controls in producing pitching moment 

Aft 

— — S exhibited a normal, trend of reduced effectiveness with increased 
A 6 

Mach number. 


Langley Aeronautical Laboratory^ 

National. Advisory Committee for Aeronautics^ 
Langley Field, Va., November 2, 1955 » 



Clarence A, Brown, Jr. 
Aeronautical Research Scientist 
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TABLE I 


The physical characteristics determined hy preflight measurements 
are as follows; 


Wing; 

sq^ ft 


Thickness/Chord at body jimcture 
Wing span^ ft o , 
Aspect ratio^ exposed 


2.839 

1,650 

0,012 

1.750 

0,854 


Canard control surfaces; 

S 0 ^ sq ft 
c ^ ft 

Thickness/Chord at body junct'ure 
Control-surface span^ ft , . . , 
Aspect ratio^ exposed 


, 0.358 
. 0,651 

, 0.033 

, 1,260 
, 1.980 


General; 

Boby dxsmeter ^ xn. ...e. 

Exneness ratxo 

A^ body cross-sectional area^ sq ft . 

Weighty lb (model sustainer loaded) 

Weighty lb (model sustainer empty) 

Moment of inertia; 

p 

Modelj, sustainer empty, ly, slug-ft 

Model, sustainer empty, slug-ft^ 
Center-of~gravity location, model sustainer empty, 
inches from nose 

Center-of-gravity location, model sustainer loaded, 
inches from nose 

Ratio of span of control surfaces to span of wings » » 


. 5 • 000 
. 22,95 

. 0.156 

158.25 

121.25 

. 3T»52 
. 0,215 

, 50.50 

, 58.00 
L 0.72 
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Figure 1.- Sketch of model 



Side view 


All dimensions are in inches 





»30 


Section A-A 



Section A-A 

Mode! wing surface 

Figure 2,- Sketch of control surface and wing surface for model tested. 

All dimensions are in inches. 
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Model rotated 


Figure 3.- Photographs of model tested 








vjr^ 


Figure 4.- Photograph of model and booster prior to laimching 
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(a) 5w-0.3°. 

Figure 6.- Variation of lift coefficient with angle of attack 
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Figiore 6.- Continued 






(c) 6 » -0,5°. 

Figure 6»- Continued, 
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Figure 
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& « +5.0'-’. 
6.- Continued 









(e) 6*«-0.3°. 

Figure 6.- Continued 
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(f) 5 «s +5-0°. 

Figure 6.- Concluded, 


CONFIDENTIAL 



CCIJFIDENTIAL 



Figure 7 ,- Variation of average lift-curve slope with Iyfe.ch number. 
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Figure 8.- Variation of the exponential damping constant h with Mach number. 
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Figure 9 .- Variation of the aerodynamic damping-in-pitch derivative 

with Nfe.ch nimiber. 
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Fig\ire 10.- Variation of period of oscillation with Bfech niuriber. 



Figure 11.- Variation of the static stability derivative Cm with Nfech 

number. 
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Figure 12. - Continued 
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Figure li|.- Variation of trim angle of attack with Wfech number. 



Figure 15*. - Variation of the trim angle of attack produced by a unit 
control deflection with tfech number. 
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Figure l6,- Variation of the lift per unit control deflection with Mach 

nuinher. 



Figiire I 7 . - Pitching effectiveness of canard control surfaces 
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Figure l8.- Variation of the drag coefficient with Mach number of the 
models of references ^ and B and the model tested. 
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